Flory-Huggins theory for the solubility of heterogeneously-modified polymers 
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Many water soluble polymers are chemically modified versions of insoluble base materials such as 
cellulose. A Flory-Huggins model is solved to determine the effects of heterogeneity in modification 
on the solubility of such polymers. It is found that heterogeneity leads to decreased solubility, with 
the effect increasing with increasing blockiness. In the limit of extreme blockiness, the nature of the 
phase coexistence crosses over to a polymer-polymer demixing transition. Some consequences are 
discussed for the synthesis of partially modified polymers, and the experimental characterisation of 
such systems. 



Many water-soluble polymers are made by chemically 
modifying insoluble base materials such as starches and 
gums, for example a wide class of water-soluble polymers 
are obtained from cellulose [l|, l2| . It is often possible to 
vary the degree of modification of the base polymer to 
obtain water soluble polymers with, in principle, con- 
tinuously variable properties. A basic characteristic of 
these polymers is their solubility, but given the essen- 
tially stochastic nature of the chemical modification step, 
what is the effect of heterogeneity in modification on the 
solubility of the resulting materials? 

In the present paper, this question is approached from 
a theoretical point of view by setting up a Flory-Huggins 
model for the phase behaviour of a polymer-solvent mix- 
ture [3J, where the polymers have a random degree of 
modification. In this approach, the issue of solubility 
is translated into the problem of determing the phase 
coexistence between a dissolved aqueous phase and an 
undissolved (water-poor) phase. The solubility is then 
formally given by the polymer concentration in the aque- 
ous phase. Determination of the full phase behaviour for 
a multicomponent Flory-Huggins theory is an onerous 
task though, and a simpler approach is to examine the 
spinodal stability of the system, which can be taken to 
be representative of the full phase behaviour. This is the 
approach taken in the present paper. It is arguably more 
insightful than a full calculation of the phase behaviour 
since closed-form analytic expressions can be obtained for 
the spinodal stability limit. The approach taken is sim- 
ilar to models for the phase behaviour of random block 
copolymer melts which have been developed in the past 
0; IE HI- There has been rather little work though on 
random copolymers which also include a solvent, apart 
from a brief example described by Sollich et al 0] ■ 

In the present model, it is supposed that the system 
comprises a large number of species of polymers i with 
differing degrees of modification < 04 < 1 and con- 
centrations pi. For simplicity, length polydispersity is 
neglected, and all the polymers are assumed to have the 
same number N of segments. The system is then de- 
scribed by the following (mean field) Flory-Huggins free 
energy density, 

/ = E I Paogp J + (l-0)log(l-0)+ X (0-r y )(l-0), (1) 

where 4> is the total polymer segment concentration and 



r\ is the concentration of chemically modified segments, 
given respectively by = Nj^iPi and 77 = N"£ ii p i a i . 
The first term in Eq. |1| is the ideal free energy of mix- 
ing. The second term is the usual Flory-Huggins configu- 
rational chain entropy. The third term is the free energy 
cost of the unmodified polymer segments at a concentra- 
tion <f> — 77 coming into contact with solvent (water) at a 
concentration 1 — 0. Typically one expects x > 1/2 for 
this interaction, to represent the repulsion between un- 
modified segments and water which leads to phase sepa- 
ration of unmodified polymers. To keep the model sim- 
ple, this is the only x-parameter that is retained in the 
problem. 

Eq. |T]) has the structure of a moment free energy, 
since the excess free energy, comprising the second and 
third terms, only depends on <f> and rj which are mo- 
ment densities. Such a system can be analysed using the 
methods developed by Sollich and coworkers @, H, l9l.ll0|| . 
In particular, Ref. 10] describes how the spinodal sta- 
bility conditions for systems with an excess free energy 
can be expressed in terms of moment densities, gen- 
eralising various truncation theorems obtained by ear- 
lier workers [TJ E3] • I now summarise the relevant re- 
sults, translated into terms suitable for the present prob- 
lem. Let us consider such a system with a free energy 
/ = J2i Pi log Pi + f^ift 1 ' ...(f)^), where the excess 
free energy depends on moment densities of the form 
0(r) _ ^2 i p i w^ %> (t = l...n), with the w^^ being 
species-dependent weights. The fundamental idea is that 
the moment densities can be treated as effective species 
concentrations. In particular, it can be proved that spin- 
odal stability corresponds to the positive-definiteness of 
the matrix M of second partial derivatives of the free en- 
ergy with respect to the moment densities. In Ref. [To| 
it is shown that M = M id + M ox where (M^ 1 )^ = 

J2iPi^i r) wi s) and (M cx ) rs = <9 2 / (cx) /<M r) cV (s) - The 
limit of spinodal stability is given by det M = 0. This 
condition usually corresponds to the vanishing of a single 
eigenvalue of M, with an eigenvector A(f>^ that satisfies 
£ a (M) rs A</>00 = 0. It is shown in Ref. [nj that the 
spinodal instability direction in the space of species con- 
centrations is given by Ap t = J2 rs Pi w i ( M id)rsA0 (s) . 

For the present problem, there are two moment densi- 
ties 4> and 77, defined respectively with = N (a con- 
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stant) and w\ = Noti (the number of modified groups 
on the ith species). Application of the above theory to 
Eq. (PP> leads to 



and 



M" 1 = N 2 




(2) 



(3) 



X 

After some algebra the condition det M = reduces to 

-L + -1^ - 2 X (1 - <«)) - X 2 iV 0«a 2 ) - (a)) = 0, (4) 
where 

( a ) = Y,iPi a i/Y,iPh (" 2 ) = E t A a ?/E 4 A- ( 5 ) 

I emphasise that, despite being remarkably simple, 
Eq. (jU is exact. 

One already reaches a significant conclusion from this. 
The first three terms in Eq. (U) are what one would ex- 
pect from standard Flory-Huggins theory Q , with an ef- 
fective x-parameter given by the product of the original 
X-parameter and the fraction 1 — (a) of unmodified seg- 
ments. These terms therefore take account of the mean 
degree of modification. The final term in Eq. (j4|) is a 
correction due to the heterogeneity. Since the variance 
(a 2 ) — (a) 2 is positive, this term is always negative. The 
effect is that heterogeneity in modification reduces the 
solubility, over and above what would be expected from 
the mean degree of modification. 

To make further progress, it is convenient to specify a 
model for the distribution of the o^. In particular, such a 
model can be used to examine the effect of blockiness in 
modification which is expected to play an important role. 
In previous work on random block copolymers 0, Q , a 
Markov model was used to characterise the correlations 
between different kinds of segments. Whilst such a model 
may be appropriate for the stochastic nature of the syn- 
thetic route for such random block copolymers, as dis- 
cussed below it is probably not appropriate in the present 
case. I therefore consider instead a very simple model for 
the heterogeneity in which the modified segments occur 
in blocks of size M, where I < M < N. In this model, 
it is supposed that each block has an equal probability p 
of being modified, and there are no further correlations. 
Then, for any particular species, on — (1/-ZV) J^flf 1 
where j labels the blocks, and is zero or one with prob- 
ability 1 — p and p respectively. Thus the a, are drawn 
from scaled binomial distribution, with 



(a)=p, (a 2 )-(a) = (M/N)p(l-p). 



(6) 



Eq. (ID) becomes 

I 1 

iV^ + 1 - " 



This is a quadratic equation for x an( i the appropriate 
root is 
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(8) 



2 x (l-p)^ 2 l%(l-p)=0. (7) 



I now examine the consequences of this result. 

The formal limit M — > corresponds to a vanishing 
variance and a completely uniform distribution of mod- 
ified segments, as though each monomer has undergone 
an identical fractional modification by a fraction p, rather 
than being modified or not with probability p and 1 — p. 
As noted already above, this limit corresponds to sim- 
ple Flory-Huggins theory with an effective ^-parameter 
equal to x(l —p)- For large N, this indicates the absence 
of phase separation for x(l —p) < 1/2 or p > I — l/(2x). 

Now let us consider Eq. {SJ for block size M = 1. 
In this case, individual segments are modified randomly 
with no correlations. For M = 1 and large N in Eq. f8j, 
there are two behaviours depending on the value of p. 
For p < 4/5, there is an absence of phase separation 
for x(l — p) < 1/2, just as for the M — > limit. For 
4/5 < p < 1, the behaviour is more complicated. To be 
precise, the location of the minimum value of the x(</0 
spinodal shifts from </> m i n ~ N^ 1 / 2 for p < 4/5 to a non- 
vanishing < m in < 1 for p > 4/5 (it is the examination 
of Eq. ([8]) in the limit cp ~ TV -1 / 2 that gives the cross over 
point p — 4/5). The change in behaviour can be seen for 
the M — 1 curves (dashed lines) in Fig. [T] and is shown 
explicitly in the upper plot of Fig. [2] 

Let us next consider the limit of extreme blockiness 
M = N. This limit is strikingly different from the M = 1 
case. For large N and p > 0, one can show that there is an 
absence of phase separation only for x\/ Np(\ — p) < 2. 
In the large N limit, this inequality is always violated, 
indicating that the system always has a tendency to un- 
dergo phase separation in the limit of extreme blockiness. 
Since the unmodified polymer system itself only phase 
separates for x > 1/2, this suggests that the phase sep- 
aration has the nature of a polymer-polymer demixing 
transition rather than a solvent-driven phase separation. 
This insight is confirmed by analysis of the spinodal in- 
stability direction below. 

For large N and general M in Eq. ©, one would ex- 
pect that the above two cases represent the two classes of 
behaviour. In the first case M <€. N and the behaviour 
is similar to the M = 1 limit where individual segments 
are randomly modified. In the second case, M oc N 
and the behaviour is similar to the M — N limit of ex- 
treme blockiness. Fig. Q] shows typical spinodal curves 
calculated from Eq. f8} for various values of p and M. 
The location of the minimum (<?!> m in, Xmin) of the spin- 
odal curves can be numerically determined, and Fig. [5] 
shows how this depends on p. 

The results show firstly that for M N, increasing p 
leads to increasing solubility as the value of x required to 
reach the spinodal instability is increased. Moreover, a 
decrease in solubility between a uniform model (M — > 0) 
with no heterogeneity, and a model with fine-grained 
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FIG. 1: Spinodal curves calculated from Eq. (0 for polymers 
of length N = 10 3 , for three values of the mean degree of 
modification p, and for block sizes M — » (uniform limit, 
solid line), M = 1 (dashed line), M = 10 (dash-dot line), 
M = 100 (dash-dot-dot line) and M = 10 3 (dash-dash-dot 
line). The system is spinodally unstable above the indicated 
curves. Note the change in shape of the M = 1 curves: for 
p — 0.1 and 0.5 the minimum is at <f> — > 0, whereas for p = 0.9 
the minimum is at <f) w 0.25. 



blockiness (M = 1), is apparent. The major effect arises 
as M — > N though, where the tendency for phase sepa- 
ration is greatly enhanced. 

The above analysis is augmented considering the spin- 
odal instability direction associated with the spinodal 
stability limit which can provide a useful mechanistic in- 
sight. As explained above, the spinodal instability direc- 
tion is characterised by the eigenvector that corresponds 
to the vanishing eigenvalue responsible for the vanish- 
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FIG. 2: The location of the numerically determined minimum 
of the spinodal curves from Eq. ([8| is plotted as a function 
of p, for polymers of length iV = 10 3 and block sizes M — > 
(uniform limit, solid line), M = 1 (dashed line), M = 10 
(dash-dot line), M = 100 (dash-dot-dot line) and M = 10 3 
(dash-dash-dot line). For M = 1 (dashed line) the upper 
plot shows clearly that m i n » N~ 1/l2 w 0.03 only holds for 
p< 4/5 = 0.8. 



ing spinodal determinant. For the present problem, from 
Eqs. ([21)- 0, one finds the instability direction is charac- 
terised by 

Ar 1 /Aq>=(a)- X N0({a 2 )-(a) 2 ) (9) 

The corresponding spinodal instability direction in the 
space of species concentrations is 

Apj _ {a 2 )Acj) - (a)Arj + a^Ar] - (a) Acj>) 
~ 0((a 2 )-(a) 2 ) 

where the second line follows by inserting the result for 
the ratio Arj/A(f). These results should be evaluated on 
the spinodal. They are all exact, for an arbitrary distri- 
bution of OLi. 

For the instability direction to lie along a pure dilution 
line, one should have Api/pi independent of species i. 
One can conclude that this only happens if Ary/A</> = 
(a), in other words if the variance (a 2 ) — (a) 2 vanishes. 
In such a case, the phase transition is purely associative, 
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or solvent-driven, meaning that the compositions of the 
coexisting phases remain the same (Api/Aip — Pi/4>)- 

If one specialises to the model of blockiness described 
above by inserting the value of \ corresponding to the 
spinodal stability limit, the instability direction becomes 
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This confirms that the spinodal instability lies along a 
dilution line (Ar/fAcf) = (a) = p) only in the limit 
M — > which formally corresponds to a vanishing vari- 
ance. For M = 1 (and M -C N in general) the phase 
transition has a mixed character. The interesting case 
occurs when M = N (or M oc N in general) for which 
Ar]/A(j) ~ (-)N 1 / 2 in the limit of large N. One can write 
this as A0/A?7 — > as N — » oo. This shows that the 
phase transition tends towards being purely segregative, 
meaning that the overall polymer concentration in coex- 
isting phases remains the same (A<p = 0). This confirms 
the suggestion above, that in the limit of extreme block- 
iness, the system tends towards a segregative polymer- 
polymer demixing transition. 

Let us now try to draw some conclusions. The main ef- 
fect of randomness is to reduce the solubility of partially- 
modified polymers beyond what would be expected from 
the mean degree of modification. The extent to which 
this occurs depends on the blockiness in substitution. 
For fine-grained blockiness, the phase behaviour is ex- 
pected to be similar to a system for which there is no ran- 
domness, albeit with a somewhat reduced solubility. For 
coarse-grained blockiness, where the block size is com- 
parable to the polymer length, the nature of the phase 
transition changes to a polymer-polymer demixing tran- 
sition. In this situation, one expects that the modified 
polymers (being almost fully modified) will partition into 
the aqueous phase, leaving the unmodified polymers be- 
hind. 

The reason for considering the two extreme kinds of 
blockiness is now clearer: namely one can envisage two 
different mechanisms of chemical modification (this is the 
reason why a Markov model for the distribution of mod- 
ified segments has not been used). Fine-grained block- 
iness would arise if monomers are equally accessible to 



the modifying agent, irrespective of their surroudings. If 
this cannot be achieved in a one-step process (for the 
reason described below) it could perhaps be achieved in 
a two-step process, by fully modifying the polymers then 
removing a random fraction of the derivative groups. Ex- 
treme blockiness on the scale of the polymer chain itself 
would arise if the modifying agent was present only in the 
aqueous phase, and as such only able to access polymer 
which had already been solubilised. This would lead to a 
mixture of polymers which were cither fully modified, or 
remained unmodified and insoluble. The process of mod- 
ification of insoluble polymers could still be initiated be- 
cause the modifying agent is able to access the tiny pro- 
portion of the insoluble polymer segments which lie at 
the interface between the insoluble and aqueous phases. 
Experimentally, confirmation of the scenario of extreme 
blockiness would be given by measuring the mean degree 
of modification for the dissolved polymers. One should 
find that this is much in excess of the apparent mean 
degree of modification. 

In the calculation, the major effect arises from inter- 
chain rather than intra-chain heterogeneities. The model 
is not sophisticated enough to take account of the solu- 
tion structures such as micelles or mesophases that could 
form for blocky polymers with block sizes M 3> 1 but 
still M < N (for example, diblock copolymers). Such 
polymers would be expected to have greater solubilities 
than would be predicted from the Flory-Huggins theory 
since the hydrophobic groups can be buried in micelles 
or other solution structures. The present theory could be 
extended to discuss these inhomogeneous situations using 
a Landau approach developed for random block copoly- 
mers [4j, |5|, |6| . For the mechanistic routes discussed above 
though, it is difficult to envisage that polymers with in- 
termediate block sizes could arise very easily. I therefore 
expect that the general conclusions will remain. 

Finally I note that in principle the above model for 
the phase behaviour could be combined with a model for 
the chemical modification reaction, to obtain a theory for 
reaction-induced solubility. However, one needs to take 
great care to capture the kinetics correctly [l3| . 

I thank Nigel Clarke for a critical reading of the 
manuscript. 
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